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parameter minima. Introduced initially for strong first-order phase 
transitions in q-state Potts models with "large enough" q, the new 
techniques prove to be surprisingly accurate for a q value as small as 
3. On the basis of the high-precision Monte Carlo data of Alves et 
al. [Phys. Rev. B43 (1991) 5846], this leads to a refined estimate of 
Pt = 0.550 565(10) for the infinite- volume transition point. 
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1 Introduction 



The three-dimensional (3D) 3-state Potts ferromagnet serves as an important 
model in both condensed matter as well as high-energy physics [1]. Experi- 
mental realizations are structural phase transitions in some crystals, and the- 
oretically this model has attracted much interest as a simple effective model 
of finite-temperature pure-gauge QCD. Consequently it has been studied in 
the past few years by many authors using quite a variety of different tech- 
niques [2-23]. By analyzing Monte Carlo (MC) simulations with the help of 
standard finite-size scaling (FSS) methods [24-27] the characteristic param- 
eters of the phase transition (transition temperature, latent heat, etc. ) have 
been estimated with varying accuracy. As a result there is by now general 
consensus that this model undergoes on a simple cubic lattice a weak first- 
order phase transition from a three-fold degenerate ordered low-temperature 
phase to a disordered phase at high temperatures. 

Most high-precision MC results are based on the FSS of the specific- 
heat maxima, the Binder-parameter minima, or the partition function zeros. 
For some of these observables the pseudo-transition points, Pt{L), exhibit a 
non-monotonic FSS behaviour with a peculiar dip around linear lattice sizes 
L = 30 for, e.g., the specific- heat data. This clearly indicates that for the 
standard observables the asymptotic FSS region is not reached until L > 30. 
For large systems with periodic boundary conditions, as a consequence of 
phase coexistence at a first-order phase transition, the FSS behaviour of these 
observables is governed by asymptotic series expansions in inverse powers of 
the volume V = [24-27]. In addit ion there are also further corrections 
that decrease exponentially with the system size [28, 29]. They originate from 
finite-size effects oc exp(— L/Lq) in the pure phases, where Lq is of the order 
of the finite (pure phase) correlation lengths ^o,d in the ordered and disordered 
phase, and from contributions oc cxp[—2aodL^~^) of the two-phase region, 
where aod is the (reduced) interface tension between the two phases. For 
3D systems the latter correction decays asymptotically much faster and, in 
general, both types of exponential corrections are always weaker than the 
power-law terms for large L. For relatively small system sizes with L of the 
order of ^o^^ or o'od^^^^^\ however, it is a priori not clear which type of scaling 
behaviour will dominate. A possible explanation for the observed dip in the 
scaling behaviour of the 3D 3-state Potts model is therefore that for L < 30 
the exponential terms dominate and then rapidly die out, giving way to an 
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inverse volume dependence for L > 30. 

Consequently, in order to extract the infinite- volume transition point, (3t, 
from a pure power-law ansatz in the inverse volume, in previous studies only 
the data from lattices with L > 30 have been taken into account. Of course, 
this is an expensive task because the computer time to generate accurate 
enough data points for such lattice sizes becomes exceedingly large, even 
with refined MC techniques for first-order phase transitions such as multi- 
canonical [30] or multibondic [31] algorithms which reduce autocorrelation 
times significantly. 

Recently a new set of observables was proposed which do not show any 
power-law terms in their FSS behaviour [28, 29]. It is thus governed only by 
exponentially small terms. For lattices with periodic boundary conditions 
and large enough g, this could be proven exactly [25]. In Refs. [28, 29] the 
new set of observables was first used to study the first-order phase transitions 
of the two-dimensional g-state Potts model with g = 5, 8, and 10. It was 
found that the new methods yield surprisingly accurate results already for 
very small lattice sizes and even for the extremely weak first-order transition 
of the 2D 5-state model. Since this might be a fortuitous accident of the 2D 
model we found it worthwhile to explore the accuracy of the new methods in 
the 3D case as well. The purpose of this paper is thus to test the validity of 
the new observables for three-dimensional lattices and for a q value as small 
as 3. We shall see that, due to the absence of power-law corrections, the dip 
in the new pseudo-transition points PtiL) indeed disappears, and a good fit 
to extract the infinite-volume transition point Pt = A(oo) becomes feasible 
without the need of extremely large lattices. 

The rest of the paper is organized as follows. In Sec. 2 we recall the 
definition of the model and briefiy summarize previous results, and in Sec. 3 
we recapitulate the definition of the new observables. The results of our 
analysis are presented in Sec. 4, and in Sec. 5 we conclude with a brief 
discussion. 

2 Model and previous results 

We use the standard definition of the 3-state Potts model [1], 

Z=^e-^^ = ^7V(£;)e-^^, E ^ -^,5,^,,^, <7, = 1,2,3, (1) 

{Ti} E (ij) 
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where = J/ksT is the inverse temperature in natural units, N{E) is the 
number of configurations with energy E, {ij) indicates that the sum runs over 
nearest-neighbour pairs, and Sfji,a- is the Kronecker symboL The energy per 
site is denoted by e = E/V, where V = is the volume of the system. The 
data we use to test our observables was obtained in a previous MC simulation 
by Alves et al. (ABV) [12] who employed heat-bath updating for lattices of 
size L = 10, 12, 14, 18, 22, 24, 30, and 36, with periodic boundary conditions. 
For all the lattice sizes, the simulations were performed at Pmc — 0.55059, 
with 50 X 10^ energy measurements for L — 10, 12, and 14, and 20 x 10^ 
measurements for the larger lattice sizes. 

By means of standard reweighting procedures [32], the data was used 
to compute the maximum Cmax = C'(/5cmax) of the specific heat, C(/5) = 
/9^y((e^) — (e)^), and the Fisher zeros of the partition function. Finite-size 
scaling was finally applied to both the specific-heat maxima and the parti- 
tion function zeros in order to extract the transition point and the FSS expo- 
nents. Due to the above described problem with the non-monotonic scaling 
behaviour of /3cmax(-^)) ABV combined their data with those of Ref. [11] on 
larger lattices and extracted their final estimate^ of 

/3t = 0.550 523 ±0.000 010 (ABV) 

from a power-law fit, l3c\^^^ = I0t + c/V , using only the largest lattices of size 
L = 30, 36, 42, and 48. The latent heat, 

Ae = e<i- Co = 0.160 62 ±0.000 52 (ABV), 

was obtained from the FSS of the specific-heat maxima on smaller lattices 

of size L = 22, 24, 30, and 36, since the scaling behaviour of Cmax turned 
out to be better behaved than that of the maxima locations. An overview of 
other previous estimates of Pt is given in Table 1. 



3 The new observables 

The definition of the new observables is based on the observation made in 
Ref. [25] that for a g-state Potts model on a lattice with periodic boundary 

^The notation of ABV follows the high-energy physics conventions. The translation into 
the present notation, usually used in condensed matter physics, reads Pt = (3/2)/3(^^^''. 
Similarly the latent heats are related by Ae = 2Ae(^^^) = 2£. 
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Tabic 1: The transition point jSt of the 3-state Potts model on the simple 
cubic lattice in units k^ — ^ — J 



A 



author(s) 



0.5348 

0.559 6(54) 

0.54749(52) 

0.550 

0.55006(91) 
0.55040(19) 

0.552 

0.550 59(2) 

0.550 62(3) 

0.550545(45) 

0.550523(10) 

0.550479(61) 

0.550565(10) 



Straley (1974) [2] 

Kim and Joseph (1975) [3] 

Miyashita et al. (1979) [4] 

Blote and Swendsen (1979) [5] 

Herrmann (1979) [6] 

Knak- Jensen and Mouritsen (1979) [7] 

Ono and Ito (1982) [8] 

Wilson and Vausc (1987) [9] 

Gavai et al. (1989) [10] 

Fukugita et al. (1990) [11] 

Alves et al. (1991) [12] 

Yamagata (1993) [13] 

this work 



conditions, provided q is large enough, the partition function can be written 
as 



.m=0 



[i + o(ye-^/^°)]. 



(2) 



where Lq < cxd is a constant, L is the linear length of the lattice, and fm{f3) is 
the metastable free-energy density of the phase m, where m = is associated 
with the disordered and m = l,...,q with the q ordered phases. It can 
be defined in such a way that it is equal to the idealized infinite-volume 
free-energy density /(/3) if m is stable and strictly larger than /(/3) if m is 
unstable. At A fm are equal. 

In Refs. [28, 29] it is discussed how Eq. (2) leads to a definition of pseudo- 
transition points f^v/v' ioT finite systems, which deviate from the infinite- 
volume transition point A only by exponentially small terms. These points 
are obtained by locating the maximum of the number- of -phases observable 



N{VX,(3) = 



l/(a-l) 



(3) 
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where a = V'/V. By inserting (2) it is easy to see that N{V.V' , (3) indeed 
counts the number of phases: q at low temperatures, g + 1 at phase coex- 
istence, and 1 in the disordered high-temperature phase. Searching for the 
maximum of N(y, V', j3) as a function of j3 amounts to locating the crossing 
points of the internal energies per site for the two lattices of different size, 

e{V,(5v/v')^e{V',(5v/v'). (4) 

As a consequence of (2) it can be proven that, for large enough g, the crossing 
points (3v/v' s-re only exponentially shifted from the infinite- volume transition 
point (3t. In contrast to the pseudo-transition points Pc^^^{L), there are no 
powers of the inverse volume involved. 

Also in Refs. [28, 29], yet another definition of pseudo-transition points, 
Pw{L), without power-law corrections was proposed which are obtained by 
measuring the ratio- of -weights, 

Rw{v,p) ^ ^ ^ E PfiiE)/ E PpiE), (5) 

^^d, E<Eo E>Eo 

and solving the equation 

Rw{V,(3)\f,=p^ = q = 3. (6) 

Here, Ppi^E) stands for the probability distribution Ppi^E) = N{E)e~^^ /Z , 
with Z and N{E) defined in Eq. (1). The parameter Eq in Eq. (5) is defined 
by reweighting the probability distribution to the point (3p where the two 
peaks of Pp{E) have equal height and then taking Eq as the energy where 
Pfjp{E) has the minimum between the two peaks. This defines still another 
sequence of pseudo-transition points f3p which, however, is not expected to 
exhibit a FSS behaviour free of 1/^ powers. 



4 Results 

4.1 Number-of-phases criterion 

From Eq. (4) it is clear that the numerical determination of I3v/v' requires 
handling simultaneously two different lattices sizes. Once they are chosen, 
searching for the minimum of 

\e{V,(3v,v')-e{V\(3y,y,)\ = Q, (7) 
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Table 2: Pseudo-transition points Pv/v for pairs of lattices of size L and L' 



L\L' 


14 


18 


22 


24 


30 


36 


10 


0.552309(24) 


0.551757(15) 


0.551421(15) 


0.5513025(75) 


0.5510784(59) 


0.5509544(50) 


12 


0.551798(41) 


0.551405(21) 


0.551142(13) 


0.5510459(72) 


0.5508701(48) 


0.5507715(56) 


14 




0.551220(30) 


0.550998(17) 


0.5509160(95) 


0.5507778(53) 


0.5507039(57) 


18 






0.550808(41) 


0.550751(18) 


0.5506733(90) 


0.5506343(63) 


22 








0.550649(66) 


0.550624(15) 


0.5506043(87) 


24 










0.550618(14) 


0.5505995(83) 


30 












0.550586(18) 



using standard numerical minimizing routines combined with reweighting 
procedures, leads to the /3y/y/ collected in Table 2 for all possible pairs of 
sizes L and L' . Here and in the following statistical error bars are estimated 
by means of the jack-knife method [33] on the basis of 20 blocks. 

The question arises on how to extract an infinite- volume transition point 
from the results given in Table 2. Considering that, as discussed in Sec. 3, 
the pseudo-transition points Pv/v a-re left with only exponential deviations 
from (5ti we first tried a fit of the form 

= A + aie-"!^, (8) 

with L' > L kept fixed. It seems evident that the bigger V and V the better 
the fit. By keeping L' = 36 fixed and using L = 12, ... ,30 (i.e., the last 
column in Table 2 without the L 10 entry) we obtain an estimate of 

= 0.550 586 ± 0.000 010 (number-of-phases criterion), (9) 

with fli = 0.0027(10), bi = 0.224(34). and a goodness-of-fit parameter of 
Q = 0.99 (corresponding to a total chi-squarcd of = 0.080 with 3 degrees 
of freedom (dof)). The goodness of the fit can be visually inspected in Fig. 1. 
Strictly speaking the parameters in (8) should still depend on the (larger) 
lattice size L'. To estimate this effect we performed a similar fit with L' — 30 
fixed, using the crossing points with the smaller lattices of size L = 12, . . . , 24. 
As a result we obtained a consistent estimate of (3t = 0.550 590(20), Oi = 
0.00312(89), 6i = 0.201(28), and Q = 0.95. 
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As far as the systematic FSS corrections are concerned it would be more 
reasonable to extrapolate the /3y/y/ along the diagonal of Table 2, where 
L'jL f=i 1.1... 1.3, such that both, L and L', are sent to infinity simul- 
taneously. From the point of statistical errors, however, this procedure 
is somewhat problematic. First, the crossing points of energies for only 
slightly different lattice sizes have naturally the biggest errors since the slopes 
de{y, P) I dp and de{V', P) / dp differ only little. This is clearly reflected in Ta- 
ble 2. Second, and more difficult to take properly into account in principle, 
the data for, e.g., Pn-^ and Pi4:-^/ig-i are correlated since both involve the 
energy on the 14^ lattice. To avoid this correlation one would have to sim- 
ulate many more different lattice sizes. We nevertheless also tried this type 
of fit and, using all 6 data points along the diagonal (with L in (8) being the 
smaller of the two lattice sizes), obtained Pt — 0.550 586(15), Oi = 0.048(16), 
and hi — 0.307(26), with Q — 0.79. It is gratifying that the extrapolated val- 
ues of Pt are in perfect agreement for all three sequences of pseudo-transition 
points, even though the approach to the infinite-volume limit shown in Fig. 1 
looks quite different, in particular for the last sequence along the diagonal. 

4.2 Ratio-of-weights criterion 

From a technical point of view, the second new definition (6) of pseudo- 
transition points, Py,, is probably somewhat easier to evaluate numerically 
since it involves only one lattice size at a time. A slight complication arises, 
however, from the fact that one has first to estimate E^^ the energy cut 
separating the ordered from the disordered phase. In order to determine 
Eq, we proceeded as follows. For a given lattice size L, we reweight the 
probability distribution of the energy to the point Pp where the ordered and 
disordered peaks are of equal height, 

-Pi, max -P/3p(-E'l,max) — Ppp .max )=P2 ,max- 

(10) 

We then fix P to Pp and search for the minimum of Pgp{E) for £"s satisfying 
-Ei,max < E < -E2,max- Our Eq IS thus defined by Pmin = Ppp{Eo). Once Eq 
is determined it is kept fixed, and we reweight again over P until Eq. (6) is 
solved for P — P^,. The behaviour of RwiY, P) as a function of P is shown in 
Fig. 2. 

The Table 3 collects the Pp and P^j introduced in Eqs. (10) and (6) for L = 
10, ... , 36. Also given are the values of Pmin, employing the normalization 
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Table 3: Table of (3p, Pmin and the corresponding 2aod, together with (3, 



T 

Li 


Pp 


P 

min 




/3 

h'w 


10 


0.550 633(12) 


0.8576(26) 


0.001536(30) 


0.552 223(10) 


12 


0.550 638(20) 


0.7911(34) 


0.001627(30) 


0.5516148(72) 


14 


0.550 585(11) 


0.7310(43) 


0.001599(30) 


0.5512068(62) 


18 


0.550 572(15) 


0.5869(81) 


0.001645(43) 


0.550 878(10) 


22 


0.550 567(12) 


0.4561(87) 


0.001622(39) 


0.550 700(10) 


24 


0.550 5521(77) 


0.3984(94) 


0.001598(41) 


0.550 6586(65) 


30 


0.550 5560(65) 


0.2291(54) 


0.001638(26) 


0.550 5965(50) 


36 


0.550 5665(72) 


0.1155(66) 


0.001666(44) 


0.550 5812(59) 



-Pi, max = -P2,max = 1, aS Well aS 

2aod{L)^-HPrrAn)/L\ (11) 

which serves as a finite-volume estimator of the (reduced) interface tension 
between the ordered and disordered phase [34] to be discussed below. 

In order to extract an estimate for the infinite-volume transition point, 
we tried again an exponential fit of the form 

/3^(L)=A + «2e-'^^. (12) 

Using the Pw{L) for L = 14, . . . , 36 from Table 3, we extract 

A = 0.550 568 1 ± 0.000 005 6 (ratio-of-weights criterion), (13) 

with as = 0.00939(93), 62 = 0.1919(74), and Q = 0.41. The data for [3^ 
and f3p together with the fit (12) are shown in Fig. 3. Within error bars the 
estimates (9) and (13) are compatible with each other, but both are slightly 
higher than the previous estimate of ABV. Notice that also the estimates bi 
and 62 are consistent, indicating that Lq = l/&i,2 in (2) is about Lq 5. This 
value is roughly one half of the correlation length estimates ^ ^ 10.2, 
obtained recently in pure phase simulations of the 3-state Potts model on a 
100=^ lattice [35]. 
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In general also the pseudo-transition points f3p (where the two peaks are 
of equal height) can be used to extract (3t. As mentioned at the end of Sec. 3, 
the FSS behaviour should be qualitatively similar to /9cmax) i-®-, one expects 
an asymptotic power series 'm.l/V. A glance at Fig. 3 shows, however, that 
the 13 p for the three-dimensional 3-state Potts model are apparently almost 
constant. In fact, by simply taking the weighted average of the six vahics for 
L = 14, . . . , 36 (i.e., performing a trivial fit Pp = Pt = const.), we obtain an 
estimate of 

(3^ = 0.550 562 5 ± 0.000 003 6 (equal-peak-height locations), (14) 

with Q = 0.16, in good agreement with (13). If we omit the L = 14 value, 
the average changes only little to Pt = 0.550 560 0(38), with Q = 0.51. 

This completely unexpected result can be understood as follows: The 
decomposition (2) of the partition function translates to the probability dis- 
tribution (formally via an inverse Laplace transformation) at Pt in a Gaussian 
approximation to 



2 



(e - eo)^PtV 



2c. 




(e - egfPtV 
2cd 



= ^V^exp 

= Ho{e) + Hd{e), (15) 

where Co (e^) and Co (q) are the energy and specific heat in the ordered (dis- 
ordered) phase at Pt, and (Hd) approximates the "ordered" ("disordered") 
part of the measured histograms.^ This implies for the ratio of heights at Pt 



H 



max 



RH{Pt)^j^ = qJ-^. (16) 

Estimates of Cq and have recently been obtained in Ref. [35] from simula- 
tions of a 100^ lattice in the pure phases at P = 0.550 53, the previously best 
estimate of Pt- If we perform a slight reweighting to our present best esti- 
mate of Pt — 0.550 565 (see the Conclusions below) and insert the resulting 
numbers (co = 28.84, — 3.210) into (16) we obtain 

RniPt) = 1.001(18), (17) 



^In practice the separation of a given histogram into Ho and Hj. is of course not 
unique. Using as dividing energy the cut Eq introduced earUer, then J2e ^o{e) = Wo and 
J2e Hd{e) = Wd- The Gaussian approximation is therefore a good approximation in the 
sense that Wo/Wd = q ai (if 
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where the error emerges from the statistical uncertainties in Co and q. The 
error due to the uncertainty in Pt is of the same order; see Fig. 4, where Rh{P) 
is plotted over a range of (5 values around (5t. We thus have the surprising 
result that for the three-dimensional 3-state Potts model the first-order tran- 
sition point f3t is not only characterized by the, on quite general grounds, 
theoretically predicted ratio-of-weights R\Y{l3t) =9 = 3, but empirically also 
by an equal-peak-height condition, Rn^Pt) = 1. This is qualitatively the 
reason why (3p happens to be so close to the infinite-volume transition point 
f3f The actually measured probability distributions for L = 36 reweighted to 
(a) (3 = f3,^{L = 36) = 0.550 581 2 and (b) (3 = f3p{L = 36) = 0.550 566 5 are 
shown in Fig. 5. Notice the sensitivity to small variations in (3 of the order 
of 2 - 3 error bars of Pw or Pp. 

The result Cq/ q ~ 9 = gr^ is presumably an accidental peculiarity of this 
particular 3D 3-state Potts model. For 2D f/-statc Potts models with q > 5 
it is exactly known that Cq ~ q at Pt, such that RniPt) ~ Q, and also for 
the 3D g-state models with q = A and 5 one finds values of RniPt) ~ 1.4 
{q = 4) and RniPt) ~ 1-8 {q = 5) (using again the estimates for Co and q of 
Ref. [35]) which clearly deviate from unity. On the other hand, in (17) unity 
is hit so accurately that it is tempting to speculate that there could be some 
hidden theoretical reason for this numerical observation (possibly related to 
the symmetry of the 3-state Potts model). It would be very interesting to 
test this possibility by studying the 3-state model on other three-dimensional 
lattices (BCC, FCC, . . . ) or in higher dimensions (4D, 5D, . . . ). 

4.3 Locations of Cmax 

Next we also reconsider the scaling of the locations Pc^a.^ of the specific- 
heat maxima, showing the peculiar dip around L — 30. As discussed in 
the Introduction a possible explanation is a crossover between exponential 
corrections oc exp(— L/Lq) and the asymptotic power series in 1/V. We 
therefore tried to fit the data with an ansatz 

Pc^M^Pt + ^ + ase-"^"^. (18) 

By following the procedure of ABV and including also the data of Ref. [11], 
but discarding only the three points with L < 14, we obtained 

Pt = 0.550 552 ± 0.000 020 {C^^ locations for L < 48), (19) 
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c = -4.1(1.3), as = 0.0309(34), and 63 = 0.207(16), with Q = 0.46. The fit 
is shown as the dashed hne on the fine scale of Fig. 6. By including also the 
L — 14: data in the fit, the numbers change only slightly to Pt — 0.550 535(12), 
c = -2.86(59), 03 = 0.0356(25), and 63 = 0.2324(82), with Q = 0.41.^ 

One interesting question was if now already the data of ABV up to L = 36 
would be sufficient to get a reliable estimate of /?f . Using their five data points 
for L = 18, . . . , 36, we indeed obtained a reasonable fit with 

Pt = 0.550 569 ± 0.000 042 {C^^ locations for L < 36) , (20) 

c = -4.9(2.6), 03 = 0.0351(53), 63 = 0.207(21), and Q = 0.71. As can be 
inspected in Fig. 6 the two fits are in fact hardly distinguishable. 

As a consequence of the additional exponential term in the fit ansatz (18) 
the estimates of f3t in (19) and (20) are somewhat higher than the result of 
ABV obtained from a pure power-law ansatz. In particular the estimate (20) 
relying only on the data of ABV is now in perfect agreement with the results 
from the new methods discussed above. It appears quite natural that the 
mixing with the data of Ref. [11] perturbs this self-consistency a little bit. 

Even though the high values of Q indicate that both fits arc statistically 
self- consistent, it is a priori clear that (18) can only be an effective model of 
the true behaviour, since higher order terms in 1/V are neglected. A simple 
way to test for this is to compare the coefficient c of the 1/V term with 
the prediction of FSS theory [29, 36], c = — Ing/Ae ~ —6.5, where Ae can 
be taken from ABV or from the analysis in the next section. We see that 
within the relatively large error bars the fitted values of c are compatible 
with this expectation, but there is definitely a trend to underestimate c. We 
have checked that the pure power-law fit used by ABV, /Scmax — Pt + c/V, 
has the same problem; we can reproduce their f5t to all digits and obtain 
c = —1.77(33). This shows that the additional exponential correction helps 
at least to drive c in the right direction. A glance at Fig. 6 then shows that 
a pure power-law fit must necessarily underestimate (3f 

How about higher-order corrections in the 1/V expansion? The FSS 

■^It should be noted, however, that the ansatz (18) admits also another fit with "unrea- 
sonable" parameters (c > and as < 0) which actually has even a slightly better Q value 
of Q = 0.43. 
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expansion in 1/V up to second order reads [29, 36] 



,^c.. ^ A - III . ^ [4 - . + 0(UV% (21) 

where As = /3tAe is the entropy jump, Ae — e^ — to the latent heat, Ac = 

Cd — Co the specific-heat jump, and Cq, Co etc. are to be taken as the (pure 
phase) expectation values at j3t. Inserting f3t = 0.550 565 and again the 
estimates of Ref. [35], we obtain the explicit FSS prediction 

= 0.550 565 - 6.502 55/y + 94 804.0/1^2 + .... (22) 

In Fig. 7 the resulting first- and second-order FSS scaling curves are displayed 
as the sohd lines on the same scale as in Fig. 6. As a rough error estimate 
the dashed lines show Eq. (21) evaluated at the error bounds of (3t, i-e., at 
0.550 565 ± 0.000 010. The statistical error coming from the uncertainties in 
the pure phase expectation values is of the same order. At first sight one 
is inclined to conclude that the second-order FSS curve tends to reproduce 
the dip in the data. At its minimum, however, the difference between the 
FSS expansions to first and second order is already quite large, and one 
may suspect that the omitted terms of 0{1/V^) would change the picture 
drastically. In fact, by extending the expansion of Pcma.^ up to the fourth order 
in 1/V [37] we clearly observe the breakdown of this asymptotic expansion for 
moderate values of V. A similar observation was recently made in Ref. [38] 
for the maxima locations of C//3^ = V{{e^) — (e)^) in the two-dimensional 
10-state Potts model, using large-5 expansions to estimate the various energy 
cumulants involved. 



4.4 Latent heat and interface tension 

As proposed in Ref. [29] the rat io-of- weights method provides at (3t also a 
natural estimator for the latent heat, 

Ae{L) = -^\n(Wo/Wa)/V, (23) 

which, similar to Pw{L), should show only exponential deviations from the 
infinite-volume limit, Ae{oo) = — Co- Graphically, the VAe{L) are just 
the slopes of lni?vi/(/?) = ln{Wo/W(i) shown in Fig. 2 which can, of course, be 
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Table 4: Finite lattice estimates for the latent heat using the ratio-of-weights 
method (with Pw{oo) — 0.550 568 



L 


Ae at /3^(oo) 


Ae at l3^{L) 


10 


0.286 50(12) 


0.30404(12) 


12 


0.250 45(13) 


0.26462(14) 


14 


0.225 76(10) 


0.236 66(10) 


18 


0.195 49(23) 


0.202 97(23) 


22 


0.179 26(28) 


0.183 80(28) 


24 


0.173 92(23) 


0.17739(23) 


30 


0.165 80(30) 


0.16734(30) 


36 


0.16343(37) 


0.16433(36) 



computed in precisely this way by numerical differentiation. It is, however, 
more convenient and numerically more stable to perform the differentiation 
in (23) explicitly. Recalling the definition (5) it is easy to see that 

Ae(L) = (e), - (e)„, (24) 

where (e),^ = y J2e>Eo EPp{E) / J2e>Eo Pp{E) is the expectation value of the 
energy computed over the disordered peak only, and (e)o is defined analo- 
gously. We computed Ae(L) in two different ways. First, at the infinite- 
volume limit of (3w{L), P = 0.550 5681 ~ Pt, and second self-consistently at 
the lattice size dependent sequence oi P — (3w{L). The resulting data is given 
in Table 4 and shown in Fig. 8. 

The fits in Fig. 8 are according to the ansatz 

Ae(L) = Ae + a^e'''^^, (25) 

where Ae = Ae(oo) is the infinite-volume limit of the latent heat. When 
P = 0.550 5681 ~ Pt is held fixed, Ae(L) follows the exponential ansatz 
down to quite small lattice sizes, and by discarding only the L = 10 data 
point we obtain the infinite- volume estimate 

Ae = 0.161 21 ±0.000 27, (26) 

04 = 0.6160(80), and 64 = 0.1611(12), with Q = 0.14. When the Ae(L) are 
evaluated at Pw{L), the asymptotic scaling region is shifted to larger values 
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of L. In order to get a fit of reasonable quality, here we had to discard the 
data for L < 14. The result of the fit, 

Ae = 0.161 60 ±0.000 47, (27) 

with 04 = 0.733(50), h = 0.1597(42), and Q = 0.11, is then compatible with 
(26), and both estimates are consistent with the previous estimate of ABV 
from the scaling of the specific-hcat maxima. 

Let us finally consider the (reduced) interface tension (Tod- As can be read 
off from Table 3, the data for L = 12, ... , 36 vary so little that none of the 
more sophisticated FSS extrapolations is apphcable. Similar to the analysis 
of Pp we therefore have again simply taken a weighted average of the seven 
values with the result 

2cTod = 0.001 63 ±0.000 02, (28) 

and Q = 0.87. A similar estimate of 2aod — 0.001 57(6) was recently obtained 
in Ref. [22]. 



5 Conclusions 

In this paper we have presented a detailed study of the first-order phase 
transition in the three-dimensional 3-state Potts model on cubic lattices with 
emphasis on recently introduced quantities whose infinite-volume extrapola- 
tions are governed only by exponentially small terms. The main results are 
the following: 

(i) The phase transition is only weakly first-order. Nevertheless the new 
techniques, originally introduced for "large enough" q, i.e., strong first-order 
phase transitions, prove to be surprisingly accurate. The expected exponen- 
tial finite-size scaling of the pseudo-transition points (3v/v' and (3^ as well as 
of the "ratio-of-weights" definition of the latent heat is clearly observed. The 
results of fits of the form (3t{L) = A ± aexp{—bL) to the various sequences 
of pseudo-transition points are summarized in Table 5. Taking into account 
the small remaining systematic uncertainties by averaging over the entries 
in Table 5 we quote as our final estimate for the infinite-volume transition 
point 

/5t = 0.550 565 ±0.000 010. (29) 
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Table 5: Summary of our various estimates of (3t 



obs. 


range 


ansatz 


Pt 


Q 


Pv/v 


L' = 36, L = 12 - 30 


A + ae-^^ 


0.550 586(10) 


0.99 


Pv/v 


L' = 30, L = 12 - 24 




0.550 590(20) 


0.95 


Pv/v 


diagonal 


Pt + ae-^^ 


0.550 586(15) 


0.79 




L = 14 - 36 


Pt + ae-"^ 


0.550 5681(56) 


0.41 




L = 14 - 36 


average 


0.550 562 5(36) 


0.16 




L = 14 - 48 


A + c/V + ae-"^ 


0.550 535(12) 


0.41 


/-^Cmax 


L = 18 - 48 


Pt + c/V + ae-^^ 


0.550 552(20) 


0.46 




L = 18 - 36 


Pt + c/V + ae-''^ 


0.550 569(42) 


0.71 



(ii) The points Pp, where the two peaks of the energy probability distribu- 
tion are of equal height, show surprisingly almost no finite-size scaling. This 
in turn implies that in this particular model the transition point is not only 
characterized by the theoretically founded "3:1- weight rule", but empirically 
to a very good approximation also by an "equal- height rule" . We emphasize 
that this could be a purely accidental situation since an "equal-height rule" is 
definitely not satisfied in two-dimensional g-state Potts models for all g > 5 
and in the three-dimensional models with q = 4 and 5. 

(iii) The peculiar dip in the finite-size scaling of the locations of the 
specific-heat maxima can be explained by a fit of the form Pcm^^iL) — c/V + 
a cxp(— &L). 

(iv) The decay constant b in the exponential term exp(— 6L) is quite 
consistently found from the various fits around b ~ 0.2. This implies Lq — 
1/6 5, which is about one-half of the correlation length ~ ~ 10 
in the pure disordered and ordered phases. A similar relation was observed 
in Refs. [40, 41] for the two-dimensional 10-state Potts model, where is 
known exactly and evidence for = was obtained numerically with high 
accuracy [42]. 

(v) Similar to the discussion of the two-dimensional 10-state Potts model 
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in Rcf. [38], we find that the asymptotic expansion in \jV for /?Cmax ill" 
behaved. Only for very large system sizes higher order terms lead to an 
improvement which, however, is quite a general feature of asymptotic expan- 
sions. 

(vi) The finite-size scaling of the "ratio-of- weights" definition of the latent 
heat follows the predicted exponential behaviour. The resulting infinite- 
volume estimate, 

Ae = 0.161 4 ±0.000 3, (30) 

was found compatible with the previous estimate of ABV from the finite-size 
scaling of the specific-heat maxima. 

(vii) The interface tension o^^ between the ordered and disordered phases 
has probably not yet reached its asymptotic scaling region for lattice sizes 
up to L = 36. Experience with three-dimensional g-state Potts models with 
q — \ and 5 suggests [43] that our estimate lOod = 0.001 63(2) is a lower 
bound on the infinite- volume limit. Still it is difficult to reconcile this value 
with the estimates oiaooiPt) > 0.0048(3) and aooif^t) = 0.0040(7) in Ref. [17] 
and the quite general stability condition (TooiPt) < 2o"od [44]. Notice that at 
least in two dimensions equality seems to hold in the last relation ( "complete 
wetting"). 
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Fig. 1 

Figure 1: Finite-size scaling of the pseudo-transition points jSy/v of the 
number- of-phases criterion together with exponential fits of the form Pv/v — 
j3t + aexp(— 6L). 




Fig. 2 

Figure 2: The logarithm of the weight ratio as a function of the inverse 
temperature for lattices of size L — 10, 12, 14, 18, 22, 24, 30, and 36. 
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Fig. 3 

Figure 3: Finite-size scaling of the pseudo-transition points of the ratio- 
of-weights method together with an exponential fit of the form Pw — Pt + 
aexp(— Also shown are the points jSp where the two peaks of the energy 
probability distribution are of equal height. 
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Figure 4: Theoretical prediction (16) for the height ratio as a function of the 
inverse temperature. The vertical lines show the infinite-volume transition 
point together with the statistical error bounds. 
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Figure 5: Histogram for L — 36 reweighted to (a) = 36) 

and (b) Pp{L = 36) = 0.550 566 5. 
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Figure 6: Finite-size scaling of the specific-heat maxima locations together 
with fits of the form /Scmax = A + c/V + aexp(— 6L). 
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Figure 7: Finite-size scaling of the specific-heat maxima locations together 
with the theoretically predicted FSS power-law expansion in 1/V to first, sec- 
ond, and third order (see text). 
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Figure 8: Finite-size scaling of the latent heat derived from the ratio-of- 
weights method. The continuous lines are fits of the form Ae(L) = Ae(oo) + 
aexp(— 6L). 
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